Fermi Pockets in a d-wave Superconductor with Coexisting Loop-Current Order 
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A d-wave superconducting phase with coexisting intra-unit-cell orbital current order has the 
remarkable property that it supports finite size Fermi pockets of Bougoliubov quasiparticles. Ex- 
perimentally detectable consequences of this include a residual T-linear term in the specific heat in 
the absence of disorder and residual features in the thermal and microwave conductivity in the low 
disorder limit. 



As a theoretical proposition, it has recently been re- 
alized theoretically^^ that there exist superconducting 
phases in which there is a stable residual Fermi surface 
and correspondingly a finite density of states in the limit 
E —¥ 0. One case in which this occurs, which will be 
discussed explicitly in this paper, is the case of coexist- 
ing d-wave superconductivity and intra-unit cell orbital 
loop current order of the type proposed^, as the un- 
derlying cause of the pseudo-gap regime in the cuprate 
high temperature superconductors. Other interesting ex- 
amples are the "pair-density wave'^ and colinear FFLO 
states,^ which involve condensation of Co oper pairs with 
non-zero center of mass momentumP^El 

Two experimental developments are suggestive that 
this sort of behavior may be relevant in real materials: 

1) Riggs et aP^have noted the existence of a residual T- 
linear term in the specific heat in high quality crystals of 
the cuprate superconductor, YBa2Cu3 07_5, with a mag- 
nitude C/T ps 2mJ/mole-K 2 that appears to approach a 
non-zero value in the limit of vanishing disorder. Related 
to this is the fact that there is a finite paramagnetic con- 
ductivity at low frequencies in the superconducting state 
in the limit T — >• for samples with very low disorder^. 
If these are borne out in future experiments, it is worth 
exploring the possibility that it is due to a residual Fermi 
surface in the superc onduc ting state. 

2) A set of anomalie d 12 * 13 ! have been reported in certain 
insulating crystals with nearly triangular lattice struc- 
tures and an odd integer number of electrons per unit cell. 
These materials show no evidence of magnetic order down 
to temperatures well below the scale of the exchange in- 
teractions, and although they are insulating, they exhibit 
several Fermi-liquid-like low temperature behaviors in- 
cluding a T linear specific heat, a T independent Pauli 
susceptibility, and a T-linear thermal conductivity. It 
is an exciting possibility that this behavior reflects the 
existence of a spin-liquid phase, but there is no consen- 
sus on exactly what type of spin liquid. One possible 
set of spin-liquids which have been proposed in this con- 
text are formally analogous to superconducting phases, 
with "spinon pairing" playing the role of electron pair- 
ing. Two particular such proposals which have been con- 
sidered in this context are the spin-liquid analogue of a 
pair-density-wave state^ and a state with d-wave spinon 



pairing coexisting with a form of loop current ordei^, 
analogous to the state considered in the present paper. 

I. FERMI POCKETS IN A CURRENT LOOP 
ORDERED SUPERCONDUCTOR 

Typically in a superconducting state, the fact that the 
superconducting gap perfectly nests the Fermi surface 
suggests the existence of a gap on the entire Fermi sur- 
face, or at most the existence of gapless nodal points (in 
2D) or nodal points or nodal lines in 3D. However, sev- 
eral states are now known which admit a residual Fermi 
surface in the superconducting state. 

To be explicit, with the cuprates in mind, we consider 
a simple single orbital model on a square lattice which 
mimics the symmetry of the loop-current order derived in 
the mean-field approximation from the three-orbitals per 
unit-cell model. The non-interacting band dispersion is 
£k and there is a superconducting gap with d-wave sym- 
metry, A k = 77d(k)|Ak| (where ij d (k) = sign(& 2 - fc 2 ,) is 
the d-wave form- factor). Both et and |Ak| respect all the 
lattice symmetries, time reversal symmetry, and (for sim- 
plicity) spin rotational symmetry. The band dispersion is 
altered, ek — > £k+>/k, in the loop current state by a mean- 
field shift in the dispersion with magnitude proportional 
to the symmetry breaking order parameter. The orbital 
loop current order wc arc considering (shown schemati- 
cally in Fig. IT]) preserves spin rotational symmetry and 
is even under (x, y) — > (—y,—x), but is odd under time- 
reversal and under the reflection (x,y) — > {y,x). From 
this it follows that Jk x ,k y — ~J—k x ,-k — —Jk y ,k x an d 
thus that the quasi-particle dispersion is 

Eh, = J k ± y/el + \A k \*. (1) 

As long as the Fermi surface crosses the zone diago- 
nal, k x = iky, the pure d-wave superconducting state 
possesses gapless nodal points at k x — k y — ±q noc i e 
and k x = —k y = ±q no d e . It is important to note 
that reflection symmetry implies J{q,q) — 0, but, in 
general, J(q, —q) is non-zero, and in particular that 

J(qrwde,-qnode) = ~ J(-q n ode, Qnode) = J where J is 

how we will parameterize the magnitude of the current 
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FIG. 1. Schematic representation of the pattern of ground- 
state currents in an orbital current ordered state. 



loop order. Thus, as shown in Fig. [2j the nodes along 
the (1,1) direction are entirely unaffected by weak loop 
current order, while the nodes along the (1, —1) direction 
are shifted away from the Fermi energy. In particular, for 

k = A[q n ode + Sk X ](l, -1) + (1, 1) with A = ±1, 



Eh = A J ± 



yj (v F 8k±) 2 + (UA^fcy 



(2) 



where vp and «a are the appropriate derivatives of 
and Ak, respectively, and . . . signifies higher order terms 
in powers of 8k. Note that one node is shifted upward 
relative to the chemical potential and the other is shifted 
down by the same amount. The result is the Fermi sur- 
face shown in Fig. [2] with two nodal points and two 
elliptical Fermi pockets. The resulting density of states 
to linear order in the energy E (measured relative to the 
Fermi energy) is easily seen to be 



N(0) 



2 (J+\E\) 



(3) 



7T H VpVA 

where the first term comes from the Fermi pockets and 
the second from the two remaining nodal points. 

It is important to stress that the quasiparticles in ques- 
tion are not the quasiparticles of a Fermi liquid, but 
rather the Bogoliubov quasiparticles of a superconduc- 
tor. For instance, one can readily see that the expecta- 
tion value of the charge of a quasiparticle with crystal 
momentum k is 



<2k = e 



_£k_ 



(4) 



where the second expression is valid for states of zero en- 
ergy (i.e. on the Fermi pocket). Thus, it is easily seen 
that the quasiparticle charge is electron-like for k outside 
the underlying Fermi surface, hole-like inside the Fermi 
surface, and vanishes wherever the Fermi pocket crosses 
the underlying Fermi surface. For small magnitude of 
J, where the spectrum can be linearized in the neigh- 
borhood of the nodal point, the average of Qk along the 
Fermi pocket vanishes, but the mean squared charge is 

[Q 2 ] =e 2 [l-0(J/A)]. (5) 

Note that at least certain quantitative aspects of these 
results may be different in more complex models of or- 
bital current ordered states. For example, it has been 




FIG. 2. Fermi surface in the first Brillouin zone of a d-wave 
superconductor with coexisting intra-unit cell loop order of 
the sort discussed in the text. The Fermi surface consists of 
two elliptical pockets along the (1, 1) direction and two nodal 
points along the (I, —I) direction. 



arguerf^ that quantum fluctuations between the four 
symmetry related classical states (of which the example 
shown in Fig. 1 is one) play an unusual role in the physics 
of the state. Moreover, aspects of the electronic structure 
associated with multiple atoms per unit cell could lead to 
new features of the ordered state which cannot be cap- 
tured in the context of single-band models. 

As a point of comparison, for spin polarized electrons 
in a ferromagnet, the analysis is very similar. Here the 
quasiparticle dispersion becomes 



Eh = ha± 



d + 



(6) 



where h is the effective Zeeman field and a = ±1 is 
the spin-polarization of the quasiparticle. Now the nodal 
points for spin-up electrons arc shifted down below the 
chemical potential, and those for spin down electrons are 
shifted up, so that there are four Fermi pockets. 



II. MEASURABLE QUANTITIES 

Here we obtain theoretical expressions for some of 
the most anomalous experimentally measurable conse- 
quences of the existence of a Fermi pocket in a super- 
conducting state with coexisting d-wave superconduct- 
ing and orbital current order. We compare the results 
to those expected for a d-wave superconductor without 
coexisting loop current order, both with or without the 
effects of weak impurity scattering taken into account. 

From the expression in Eq. [3j it follows that the quasi- 
particle contribution to the specific heat at low T is 
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where the first (leading order) term comes from the pock- 
ets and the second (subleading) term from the point 
nodes and . . . signifies higher order terms in powers of 
k B T/\J\. (Here C(3) = 1.202 ... is the Reiman Zeta func- 
tion.) In the absence of orbital current order and i n the 
absence of disorder, 7 vanishes as T — > 0. Howevej 16 * 17 !, 
in the presence of disorder (at least to the extent that 
scattering can be treated in t-matrix approximation), the 
nodal density of states is broadened on the scale of the 
elastic scattering rate, %/t and consequently for T — > 0, 
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In comparing the effect of disorder with that of loop cur- 
rent order, note that ji 00 p is not only independent of r in 
the weak disorder limit, it is parametrically larger than 
Inodai by a factor 0(Jt/K). 

A well known signature of a nodal superconductor, at 
least in the strongly Type II limit (in which H c \ is neg- 
ligibly small), is the singular magnet ic-fiel d dependence 
of the low temperature specific heall 18 * 19 ! While even 
in a fully gapped superconductor there is a non-analytic 
contribution to the specific heat proportional to \H\, i.e. 
to the density of vortices, there is an even stronger field 
dependence in a nodal superconductor due to the do- 
plar shift of the quasiparticle spectrum produced by the 
superfiuid flow at long distances from the vortex core. 
Consequently, in a pure d-wave superconductor, 



lnodal{H) - Jnodal(0) 
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where (f>o = hc/e. However, in the presence of scat- 
tering, this strong field dependence is suppressed^ for 
\H\ < ^/(vft) 2 . Thus, if a residual 7 is produced by 
impurity scattering, a corresponding rounding of the \fii 
field dependence of the specific heat is expected as well. 

In a d-wave superconductor with coexisting loop cur- 
rent order, the leading order contribution to 7 is the |if |° 
term, already computed in Eq. [7J which comes from the 
pockets. However, the next to leading term comes from 
the remaining nodal points, and is thus of precisely of the 
same form as in Eq. [9j but with half the magnitude. A 
substantial residual value of 7 as both H and T tend to 
zero, coexisting with a \[H contribution to 7 which per- 
sists to asymptotically low fields would constitute strong 
evidence of loop current order (or a relative) . 

The non- vanishing density of states in the loop ordered 
state has consequences for transport properties as well. 
In general, we expect the thermal conductivity to obey 
the relation 

Knodal = (1/2)C(T) < V 2 > T = T 7 < V 2 > T, (10) 

where < v 2 > is the mean-squared quasiparticle veloc- 
ity. For the pure d-wave superconductor with weak scat- 
tering, the factor of 1/r in 7 cancels the explicit fac- 
tor of r i n the expression for k producing a "universal" 
expressio n 1 20 ! 21 1 for k, i.e. one that is independent of r. 



The parametrically larger magnitude of 7 in the loop or- 
dered state results in a parametrically larger magnitude 
of k in the clean limit of the loop ordered state 
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Note that T; oop must be calculated for the loop current 
state and not for the nodal points. 

The conductivity may be estimated by an approximate 
version of the Wiedemann-Franz law 
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apart from possible different Fermi-liquid renormaliza- 
tions and vertex corrections. Extension of this at T — > 
but finite frequencies should lead to an expression for the 
microwave conductivity 



ct(w) 



a(0) 



1 + ^Voop' 



(13) 



The integral over oj of this should then subtract from 
the superfiuid density or A^ 2 , where Al is the London 
penetration depth. 

An obvious issue is whether the existence of Fermi 
pockets in the superconducting state produces anything 
resembling the periodic oscillations as a function of 1/if, 
as would be expected in a Fermi liquid with similar dis- 
persion, i.e. with a frequency, H* = (hc/2ire)A 1 where 
A is the area of the Fermi pocket in the BZ. On the other 
hand, as discussed above, the charge changes sign on dif- 
ferent portions of the Fermi surface such that the average 
charge of the quasiparticles is close to 0. However, An- 
dreev scattering at the edges of the pockets where the 
charge is mixes the particle and hole states such that 
closed orbits in a magnetic field are in principle possi- 
ble. Evidence that such oscillations occur has been re- 
ported by Kallin et aP in the case in which the Fermi 
pockets arise in a pair-density wave. For pockets arising 
from orbital loop current order, Wang and VafekP^l have 
found oscillations of the specific heat with a frequency 
roughly proportional to A, although the field and tem- 
perature dependence appears rather more complex than 
for a Fermi-liquid. On the otherhand, Senthil et aP^, 
in a different formulation of the problem of magneto- 
oscillations in the superconducting state do not find such 
oscillations. Analytic understanding of these results is 
presently lacking, although this is clearly an important 
and interesting problem. 



III. DISCUSSION 

The existence of a well defined Bogoliubov Fermi sur- 
face in certain exotic superconducting states is a rela- 
tively new theoretical development. What is clear is that 
it likewise gives rise to complex emergent behavior at low 
temperatures, including characteristic variations of the 
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specific heat (and, presumably, other measurable quanti- 
ties) which reflect the low energy electronic structure. 

We conclude with a few comments concerning the pos- 
sible relation of our present results with experiments 
in the cuprate superconductor, YBCO. As mentioned 
in the introduction, one puzzling feature of the spe- 
cific heat in YBCO is that there is a non-zero value 
of 7(0, 0) « 2mJ/moleK 2 in the superconducting state 
of underdoped materials which does not appear to de- 
pend greatly on sample quality, and so may be intrinsic. 
This is compatible with the existence of a Fermi sur- 
face of BdG quasiparticles. If this is interpreted as being 
the result of coexisting loop order, then taking values of 
vp 2 x 10 5 to/s and v^/vp k, 1/8 deduced in experi- 
ments, one can work backward from this to obtain a value 
for J ss 60meV from the expression in Eq. [7] 

We note that the anomalous paramagnetic 
conductivity^ observed in the superconducting state 
may also be compatible with this explanation. In the 
superconducting state an integrated paramangetic con- 
ductivity S of about 2± 1 /2% of the normal state value is 
observed in YBa2Cu30g.5, which decreases to l/2±l/2% 
for YBa2Cu306.gg. S is proportional to the number of 
electrons in the pockets, i.e. S ~ {\/2ir){a 2 J 2 /Ti 2 vpVs\). 
With the J inferred from the specific heat, we would 
estimate it to be about 3% and would expect it to go to 
as loop order disappears in the phase diagram. The 
significance of the microwave conductivity experiments is 
also that the finite temperature changes of the integrated 
weight parallel the measured change of the superfluid 
density with temperature, which is linear. This puts 
additional weight to the argument that the residual mi- 



crowave conductivity and specific heat are not extrinsic 
effects due to impurities. Assuming that surviving Fermi 
pockets are, indeed, the origin of the residual specific 
heat and the uncondensed spectral weig htP, one infers 
that Ti 00 p rs 10 -11 sees. This, in turn, would imply a 
huge enhancement of the thermal conductivity over the 
"universal value" for similar tempratures and samples. 
This does not appear to have been seen, so faif22. 

There are many further experimental tests that can 
be implemented to test for the existence of near-nodal 
pockets in the superconducting state. In terms of bulk 
measurements, the most promising are systematic mea- 
surements of the field independent specific heat and para- 
magnetic conductivity at low temperature in the super- 
conducting state as a function of doping. Turning to sur- 
face probes, high resolution ARPES measurements of the 
Fermi-surface and STM studies of quasiparticle interfer- 
ence patterns could, potentially provide the most direct 
evidence for or against the existence of nodal pockets in 
the superconducting state. Note that this is an entirely 
separate issue from the issue from the origin of the "Fermi 
arcs" in the "pseudo-gap" state above T c . 
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